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Abstract 

Shell structure in the neutron-rich Ca and Ni nuclei is investigated by the spherical Hartree- 
Fock calculations with the semi-realistic NN interactions. Specific ingredients of the effective 
interaction, particularly the tensor force, often play a key role in the Z dependence of the neutron 
shell structure. Such examples are found in = 32 and = 40; = 32 becomes magic or 
submagic in ^^Ca while its magicity is broken in ^"^Ni, and = 40 is submagic (though not magic) 
in ^^Ni but not in ^"^Ca. Comments are given on the doubly magic nature of ''^Ni. We point out 
that the loose binding can lead to a submagic number A^ = 58 in ®^Ni, assisted by the weak pair 
coupling. 
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I. INTRODUCTION 



The shell structure, which is typically manifested in the magic numbers, is one of the 
fundamental concepts in the nuclear structure physics. The shell structure of nuclei is of 
importance also in astrophysics; e.g. it provides the waiting point of the s- and r-processes. 
As abundant experimental data have been obtained in unstable nuclei, it has been clarified [l| 
that the shell structure may depend on Z or more strongly than expected from most 
conventional theories. As well as the disappearance of the N = 8 and 20 magic numbers, 
the new magic numbers A^ = 16 and 32 have been indicated in neutron-rich nuclei 

a a. 

This discovery has stimulated to reexamine and refine theories with respect to the nuclear 



shell structure. The new experimental facilities j4| are expected to access heavier unstable 
nuclei in coming years. It is desirable to give predictions on the shell structure from the 
refined theories, which could be a good guidance to new experiments and will eventually be 
tested by them. 

Concerning the Z oi N dependence of the shell structure (which is sometimes called 
"shell evolution"), two mechanisms have been argued. The absent or low centrifugal barrier 
in low-i orbits may influence the shell structure near the neutron drip line 121]. Whereas 
the A^ = 8 magic number is eroded because of this mechanism, there has been no clear 
evidence for new magic numbers owing directly to the loose binding. Since the nuclear 
shell structure is formed under the average field composed of the nucleon-nucleon (NN) 
interaction, the effective A^A^ interaction may also affect the shell structure. In particular, it 
has been pointed out that the tensor force plays a significant role in the Z or N dependence 
of the shell structure sl. For full understanding of the shell structure in unstable nuclei, it 
will be necessary to take both possibilities into account. 

The mean-field (MF) theories provide us with a good tool to study the nuclear shell struc- 
ture from the nucleonic degrees of freedom. While it is yet difficult to describe structure of 
medium- to heavy-mass nuclei with the fully microscopic A^A^ interaction to good accuracy, 
the author has recently developed semi-realistic A^A^ interactions js-S], in which the Michi- 
gan 3-range Yukawa (M3Y) interaction ^ is modified so as to reproduce basic observed 
properties such as the saturation and the is splitting. The longest-range part of the central 

(C) 

channels is maintained to be the central force in the one-pion exchange potential VQp^p. 
The tensor channels in the M3Y-Paris interaction jl^ are contained in the parameter-set 



M3Y-P5' without any change. Since the significant part of the tensor force comes from the 
pions, M3Y-P5' takes well account of the leading-order effects of the chiral symmetry break- 
ing. The tensor channels are dropped in the set M3Y-P4', which is useful to investigate role 
of the tensor force. 

Shell structure of the neutron-rich Ca and Ni nuclei is an interesting topic. The = 32 



A shell model 



new magic number has been indicated by the experiments in ^^Ca 
calculation suggests that magic nature is stronger in = 34 than in = 32 [12| because of 
the tensor force, although the data on ^^Ti show no signature of the = 34 magicity jl3 |. 
Whereas A^ = 40 bef 
7, 



been given for Ca 



laves like a magic number in ^^Ni jl^, contradictory predictions have 



15| . The Z = 28 magicity has been argued in ''^Ni [l|, |5|. It could 



also be interesting whether a new magic or submagic number exists beyond A^ = 50 in the 
Ni isotopes. In this Communication we shall investigate shell structure of the neutron-rich 
Ca and Ni nuclei by applying the self-consistent Hartree-Fock (HF) calculations with the 
semi-realistic A^A^ interactions. 



II. EFFECTIVE HAMILTONIAN 



Our effective A^A^ interactions have the following form, 



(C) , (LS) , (TN) , (DD) 



= E (e^^^^TE + e^«^PT0)/(^^Hn,) L., . (S. + S,) , 

n 

,,(TN) _ / (TNE) p I .(TNO)p w(TN)/ \ ^2 q 
n 

r = (f "^^SE • [P(r.)]"'^^' + t(™)PTE ■ [P(r.)]"'^^') 5(r.-) , (1) 



where Tij = - r^, rij = \rij\, pij = (p^ - Pj)/2, Lij = Tij x p^^, Sij = 4 [3(sj ■ rij){sj ■ 
Tij) — Si ■ Sj], Tij = Tij/vij, with i and j representing the indices of nucleons, and p(r) is the 
nucleon density. Pse, -Pte, -Pso and Ppo denote the projection operators on the singlet-even, 
triplet-even, siiiglet-odd and triplet-odd two-particle states. In the M3Y-type semi-realistic 



interactions [gHsI, the Yukawa function /^^•'(r) = e '^^ ''^/pi^^'^ is employed (X = C, LS and 
TN). The density-dependent contact force 

^(DD) ig 

introduced to realize the saturation. The 
parameter-sets M3Y-P4' and P5' are presented in Ref. jsj. We note again that M3Y-P4' 



contains fopEP but with assuming v^™'' = 0, while both fQPEp and v'^™'^ of the M3Y-Paris 
interaction are untouched in M3Y-P5'. For comparison, we use the DIS parameter-set 16|] of 
the Gogny interaction, in which fl{"\r) = e~*^^" the contact form for v^^^\ and v^™^ = 
are adopted. 

It is reasonably assumed that the spherical symmetry holds in the neutron-rich Ca and 
Ni nuclei. Although the quadrupole deformation cannot always be discarded for precise 
studies, we focus on the spherical shell structure in this Communication, which is crucial 
to understand structure of these nuclei. The spherical HF calculations are implemented by 



using the Gaussian expansion method 



17H20l] and adopting the Hamiltonian H = + 



Vc — Hera., where i^Ar(= Z]jP?/2M + I]i<j%)) and -ffc.m. denote the effective nuclear 
Hamiltonian, the Coulomb interaction and the center-of-mass Hamiltonian, respectively. 
The exchange term of Vc is treated exactly. Both the one- and the two-body terms of Hem. 
are subtracted before iteration. 

It is noted that, although the DIS interaction does not contain Wopep y^™) explicitly, 
a part of their contribution is incorporated in the other channels in an effective manner. The 
same holds for M3Y-P4' that lacks v^™\ It has still been recognized j^, Q, [21I, I22I that the Z 



or N dependence of the shell structure is difficult to be described without explicit inclusion 
of fopEp and t>(™). 



III. RESULTS AND DISCUSSIONS 



For the Ca and Ni nuclei, the main correlations beyond the spherical HF solution should 
be the neutron pairing. Therefore the pair energy is a good measure for the neutron shell 
closure. The difference between the HF and the Hartree-Fock-Bogolyubov (HFB) energies 



has been presented in Refs. [3, Isl. In Ref. [15|] the neutron pairing gaps have been shown 



for the Skyrme energy density functionals SLy4 and SkM*. We further investigate the shell 
structure of the neutron-rich Ca and Ni nuclei, particularly the magic or submagic numbers 
of N, based on the spherical HF results. 




FIG. 1. Enij) of the Ca isotopes. Blue, green and red lines represent the results with the DIS, 
M3Y-P4' and P5' interactions, respectively. Dashed lines are used for positive-energy levels. 



A. Single neutron levels 



The neutron single-particle (s.p.) energies en{j) around the Fermi level are depicted in 
Fig. [1] for the Ca isotopes, and in Fig. [2] for the Ni isotopes. The HF results obtained 
from DIS, M3Y-P4' and P5' are compared with one another. To keep the figures visible, 
Enij) is shifted by a linear function of in the vertical axes, so that the Fermi energies 
do not largely deviate from the origin. The dashed lines indicate positive Enij), which may 
correspond to the single neutron resonance and is shown for reference, although the correct 
boundary condition for the resonances is not taken in the present calculations. 

n 

Because of the small difference of the symmetry energy |8|, the slope of en{j) in the DIS 
result is slightly less steep than those in the M3Y-P4' and P5' results. A notable point is 
that en{0gg/2) significantly depends on the interactions in the Ca isotopes, but not in the Ni 
isotopes. It is also noteworthy that the neutron shell structure above = 50 in the highly 
neutron-rich Ni isotopes is different from the (3 stable region. The level sequence is 1(^5/2, 
2si/2, lc^3/2 and 0(^7/2 from the lower to the higher. 
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It is helpful to view Z dependence of the neutron shell structure in order to pin down 
what gives rise to the difference between Ca and Ni. In Fig. 13 of Ref. [7|, the single neutron 
energies relative to 1^3/2, = ^n{j) — £^^(1^3/2)5 have been depicted for the = 32 

isotones as a function of Z, calculated with the interactions M3Y-P4 and P5. These Z\e„(j) 
values are relevant to the = 32 magi city. We here display Aejyj) = Snij) — £n{^P3/2) for 
j = 0/5/2 and lpi/2 calculated with M3Y-P4' and P5' in Fig. [3l in comparison with those 
with DIS. Because of the level inversion, the = 32 shell gap corresponds to z4e„(0/5/2) 
in 6°Ni, but to Z\e„(lpi/2) in ^^Ca, for the M3Y-P4' and P5' results. In the DIS case 
Aen{lpi/2) represents the shell gap both in ^°Ni and ^^Ca. Contributions of and fopEP 
to Zi^n (0/5/2) in the M3Y-P5' result are also presented. Because we are interested in the 
Z dependence which cannot be compensated by the other channels, the and ^opep 

contributions are shifted by their values at Z = 28. 

With the M3Y-P5' interaction we obtain considerable Z-dependence in Z\£:„(0/5/2). This 
is already recognized by comparing Figs. [T] and [2J This Z-dependence originates in fopEP Hi 
a„d „(™> 0. both of which act attracfve.v o„ „0/.,. as pO/,,. . occupied. v,a the mech- 
anism discussed in Refs. [5|, 121[. Note that, though the same mechanism is present also 




FIG. 3. Z\e„(0/5/2) = en(0/5/2) - en(l?'3/2) (solid lines) and Z\e„(lpi/2) = e„(lpi/2) - en(lp3/2) 
(dot-dashed lines) for the = 32 isotones. See Fig. [T] for conventions of colors. Thin red solid 
and dashed lines represent relative contributions of v^™'' and ^opep to Z\en (O/5/2) in the M3Y-P5' 
result, after shifting by their values at ^'^Ni. 



for nlpi/2, the effects are much smaller and do not lead to significant Z dependence. Not 
including these parts explicitly, the DIS interaction does not provide strong Z-dependence 
in Zie„ (0/5/2). Having fopEp but not v^™\ M3Y-P4' gives moderate Z- dependence. The 
shell gaps in ^^Ca are comparable among all the interactions. The pair energies shown 
m Refs. 0, Q confirm that ^^Ca is nearly a doubly- magic nucleus, as is consistent with 
the measured Ex{2f) value On the contrary, the experimental data in ^*^Ni show no 

enhancement of Ex{2i) 23[, suggesting meltdown of the = 32 magicity. The effects 

(c) 

of the tensor force on the mean fields, together with those of fopEp^ "well account for the 
Z-dependence of the N = 32 magicity. 

Unlike the = 32 shell gap, the difference between £n(0/5/2) and £n(lpi/2) is not remark- 
able at Z = 20 in the present calculations. As a result ^''Ca has a certain amount of the pair 
excitation, as shown in Ref. 8|]. It is emphasized that this consequence is obtained even with 



M3Y-P5' that includes reasonably strong tensor force. Thus the = 34 magicity cannot 
be concluded only from the tensor force, and influence of the other parts of the interaction 
{e.g. the central channels) on £„(0/5/2) — Sni)-'Pi/2) is important as well. 



iV = 40 



As viewed in Refs. [7|, |8|, ll5|], the magic or submagic nature of = 40 predicted by the 
MF calculations significantly depends on the input effective interactions. While the pair 
excitation is hindered both in ^^Ni and ^°Ca if we use SLy4, DIS or M3Y-P4', there is 
no signature of the A^ = 40 magicity with SkM*. All of these interactions do not contain 
the explicit tensor force. If we apply M3Y-P5' that contains realistic tensor force, the pair 
excitation is highly suppressed in ^'^Ni but not in ^"^Ca. These results are traced back to 
the s.p. energy of nOgg/2 relative to and nlpi/2. We present the Z dependence 

of AsnU) = ^n{i) - £n(0/5/2) (j = ^g<3/2 and lpi/2) for the A^ = 40 isotones in Fig. H 
The A^ = 40 shell gap is represented by Aen(S^g9/2) when Zie„(lpi/2) is negative, and by 
Aen{^gQ/2) - Aen{lpi/2) when Aen{lpi/2) is positive. 

As mentioned above, the Z\£:„(0(79/2) values obtained with DIS, M3Y-P4' and SkM* do 
not strongly depend on Z. It has experimentally been suggested that ^^Ni looks like a 
doubly magic nucleus jl^. The DIS interaction, which gives strongly suppressed pairing, 
describes Ex{2i) and B{E2) of ^^Ni to good accuracy within the quasiparticle version of the 



random-phase approximation (RPA) [26| . The sizable pair excitation with SkM* is ascribed 
to the small Aen{0gg/2). It would be difficult to reproduce the experimental data on Ex{2f) 
and B{E2) in ^^Ni with this small Aen{0g9/2)- If we use M3Y-P5' that includes realistic 
tensor force, the A^ = 40 shell gap is comparable to those with DIS and M3Y-P4' in ^^Ni, 
but Aen{0gg/2) significantly depends on Z. This Z-dependence is predominantly carried by 

I— n (C'*^ 

the tensor force, as clarified in Fig. |H Contribution of f qpep is not important in this case. 
The A^ = 40 shell gap comes minimum at Z = 20, and this leads to sizable pair excitation 



in ^°Ca jsj. As long as we rely on the shell gap in ^^Ni, it is likely that the A^ = 40 magicity 
is significantly broken in ^°Ca, since the tensor force certainly exists in the A^A^ interaction. 
This prediction will be tested by future experiments on systematics of the binding energies 
and/or of the first excited states. 

There have been arguments on the A^ = 40 magicity in ^^Ni 



25 



271, 



28|. Although the 
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FIG. 4. Aen{0gg/2) = £,1(059/2) - £,1(0/5/2) (solid lines) and zie„(lpi/2) = £„(lpi/2) - £n.(0/5/2) 
(dot-dashed lines) for the = 40 isotones. For thick lines, the same colors as in Fig. [1] are 
used. Thin black lines are the results of the SkM* interaction. Relative contribution of v^™^ to 
Z\en (0^9/2) in the M3Y-P5' result is shown by the thin red line, with shifting by its value at ^"^Ca. 



pair energy is quite small, the neutrons are still in the superfluid phase in the HFB results 
of ^^Ni, with any of DIS, M3Y-P4' and M3Y-P5'. Moreover, the experimental data show 
that the magicity is lost quickly as Z departs from 28 23j. This situation reminds us of 



the protons in ^^^Gd 



2J], and it could be more reasonable to call = 40 around ^®Ni a 



suhmagic number rather than a new magic number. 



e deformation in 



25|, by using the 



A small shell gap may induce quadrupole deformation. The quadrupol 
the iV = 40 isotones, as is known for ^°Zr, has been investigated in Ref. 
DIS interaction. Though beyond the scope of this paper, it will be of interest to study 
deformation effects with the semi-realistic interactions. 
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TABLE I. E^{2f) and B{E2;2'l 0|) in "^^Ni, predicted by the HF+RPA calculations. 







DIS 


DIM 


M3Y-P4' 


M3Y-P5' 


B{E2;2+ - 


(MeV) 
> 0+) (e^fm^) 


3.15 
84.4 


3.00 
83.4 


3.28 
87.6 


3.25 
84.4 



D. 78Ni 

In the present MF calculations, the magic number iV = 50 is maintained both for ''°Ca and 
''^Ni, irrespective of the effective interactions. The shell gap between nOgQ/2 and the upper 
levels is large enough to prevent the neutrons from being excited in the HFB calculations. 

It has been suggested [l], H] that the Z = 28 magic nature could be eroded in ''^Ni, 
because pO/5/2 comes down via the attraction from the protons occupying pO/7/2. In the HF 
calculation with M3Y-P5', such attraction is realized because the tensor force is included, 
and pO/5/2 becomes the lowest unoccupied proton orbit. However, it is not sufficient to 
violate the Z = 28 shell gap, which amounts to 5.8 MeV, via the pair excitation. Since the 
magic nature is usually linked to properties of the first excited state, we present the values 
of Ej.{2i) and B{E2) predicted by the HF+RPA calculations in Tabled As well as those of 
DIS and the semi-realistic interactions, the results of the new parameter-set of the Gogny 



interaction DIM 29| are displayed. Comparison with future experiments is desired. 



E. = 58 

References [sl, Q show that the pair correlation is greatly suppressed in ®^Ni, suggesting 
the submagic nature of = 58 near the neutron drip line. To examine the neutron shell 
structure around A^ = 58, Zi£:„(j) = en{j) — en(l'^5/2) is depicted for j = 2si/2 and 1^3/2 
in Fig. |5l Basically the interval between 2si/2 and 1^3/2 corresponds to the A^ = 58 shell 
gap. However, being mostly positive while not satisfying the correct boundary condition, 
£:„(l(i3/2) would not precisely represent resonances. Nevertheless the calculated energies of 
are useful in interpreting the current HFB and RPA results of ^^Ni, in which influence 



of the continuum is efficiently taken into account 18 



In this region we do not find remarkable interaction-dependence in the neutron shell 
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FIG. 5. Z\e„(2si/2) = e„(2si/2) - e„(ld5/2) (solid lines) and Aen{ld-^/2) = £11(1(^3/2) - en(l'^5/2) 
(dot-dashed lines) for the = 58 isotones. See Fig. [T]for conventions of colors. 

structure. As approaching the neutron drip hne {i.e. for decreasing Z), the lower-^ orbit has 
relatively lower energy because its wave function feels the weaker centrifugal repulsion and 
thereby easily extends in the coordinate space. The main correlation which may break the 

= 58 shell gap is the pair excitation out of n2si/2 to 721^3/2. The coupling between these 
two orbits via the pairing is not strong, primarily because their degeneracy (2j + 1) is small. 
Indeed, the coupling matrix element J = 0\vij\{n2si/2y J = 0) is ~ 0.3 MeV if 

evaluated by the M3Y-P5' interaction, appreciably smaller than 2[£„(l(i3/2) — en{2si/2)] ~ 
3 MeV. Thus, the loose binding, assisted by the weak coupling, leads to the = 58 submagic 
nature in *^Ni, although the pair excitation remains within the HFB regime. 

Because the weak coupling plays a certain role, the submagic number = 58 at ^^Ni 
does not imply high Ex{2f). On the other hand, it could be manifested by suppressed 
2f — > Oi"). The 2f state is not easily handled in numerical calculations, because it is 
located just above the neutron threshold. Though the boundary condition should be treated 
more carefully to get precise values, we note for reference that the HF+RPA calculations 
using the Gaussian expansion method [2^ give Ex{2i) ~ 1 MeV and B{E2;2i — )■ Oi) = 
10 -20eW. 
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IV. SUMMARY 



We have investigated shell structure of the neutron-rich Ca and Ni nuclei by the spherical 
Hartree-Fock calculations mainly with the semi-realistic A^A^ interaction M3Y-P5'. In Z 
dependence of the neutron magic or submagic numbers, specific ingredients of the effective 
interaction, particularly the tensor force, could play a crucial role. The magic nature of 

= 32 around ^^Ca and the non-magic nature around ^°Ni can be accounted for by the 
tensor force as well as by the central part of the one-pion exchange potential; i.e. the 
leading order effects of the chiral symmetry breaking. On the other hand, the present mean- 
field study does not support the = 34 magic number. The tensor force gives rise to Z 
dependence of the shell structure around A^ = 40. Whereas the submagic nature of A^ = 40 
in ®^Ni has been observed and is described by many effective interactions except SkM*, 
the submagic nature is likely destroyed in ^°Ca because of the Z dependence in the shell 
structure produced by the tensor force. 

Although it has been pointed out that the loose binding could lead to new magic numbers 
in drip-line nuclei, no clear evidence has been found so far. We point out that A^ = 58 will 
be submagic in ^^Ni, owing to the lower centrifugal barrier in the lower I orbits, together 
with the weak pair coupling. This submagic nature may be connected to small B{E2), but 
not to high E^{2t). 

In the present work we have constrained ourselves to the spherical MF calculations, which 
are useful to understand variation of the structure in the Ca to Ni nuclei in a simple manner. 
Future plan includes extension of the calculations by taking the possibility of deformation 
into consideration |l9l| . as has been done with the phenomenological Skyrme or Gogny energy 
density functionals IsQ]. 
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